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Validity of local equilibrium has been questioned for non-equilibrium systems which are char¬ 
acterized by delayed response. In particular, for systems with non-zero thermodynamic inertia, 
the assumption of local equilibrium leads to negative values of the entropy production, which is in 
contradiction with the second law of thermodynamics. In this paper we address this question by 
suggesting a variational formulation of irreversible evolution of a system with non-zero thermody¬ 
namic inertia. We introduce the Lagrangian, which depends on the properties of the normal and the 
so-called ’’mirror-image” systems. We show that the standard evolution equations, in particular the 
Maxwell-Cattaneo-Vernotte equation, can be derived from the variational procedure without going 
beyond the assumption of local equilibrium. We also argue, that the second law of thermodynam¬ 
ics should be understood as a consequence of the variational procedure and the property of local 
equilibrium. For systems with instantaneous response this leads to the standard requirement of the 
local instantaneous entropy production being always positive. However, if a system is characterized 
by delayed response, the formulation of the second law of thermodynamics should be altered. In 
particular, the quantity, which is always positive, is not the instantaneous entropy production, but 
the entropy production averaged over the period of the heat wave. 


I. INTRODUCTION 

In classical irreversible thermodynamics (CIT), heat 
propagation is described by Fourier’s law, which states 
that the heat flux is proportional to the temperature gra¬ 
dient, J = — XVT. This relation is in agreement with the 
second law of thermodynamics. Indeed, in case of heat 
conduction, the entropy production is proportional to the 
product of the heat flux and minus the gradient of the 
temperature pj. Thus, the entropy production is always 
positive (zero in equilibrium), which is one of the state¬ 
ments of the second law of thermodynamics. 

However, Fourier’s law of heat conduction is not the 
only known constitutive relation Hi- In particular, 
some systems can exhibit delayed response 0 or even 
wave behavior 0. A common extension of Fourier’s law 
of heat conduction is the so-called Maxwell-Cattaneo- 
Vernotte (MCV) equation 0-0 

r j + J = —AVT (1) 

where r is a parameter with the dimensionality of time, 
and dot above a variable indicates the partial time deriva¬ 
tive. The modification with respect to the standard 
Fourier equation is the term r J. It can be viewed as 
the first term in the Taylor expansion of the flux J(f + r) 
around J(t), so the entire left hand side of Eq. JTJ rep¬ 
resents the delayed response of the system at the time 
t+r to the perturbation at the time t. Delayed response 
implies that the system has some sort of thermodynamic 
inertia, and the parameter r is a mere of this inertia. 

Existence of an additional term proportional to the 
temporal derivative in the evolution equation changes the 
nature of the temperature evolution. The system is no 
longer purely relaxing. In particular, if r is large enough, 
the second term on the left hand side of Eq. CD may be 
neglected and we will obtain a wave equation. For not 


very large but non-zero values of r the system behaves 
similarly to a damped harmonic oscillator Q. In partic¬ 
ular, it relaxes towards equilibrium, and, depending on 
the values of the parameters r and A, may also exhibit 
oscillating behavior. Another important consequence of 
existence of the first term in Eq. HD is the restriction 
on the speed of heat propagation. In a system described 
by Fourier’s law, a thermal perturbation propagates with 
infinite speed, which is one of the major critics of CIT. 
However, in the system described by Eq. HD this speed 
is finite and is proportional to 1/^/r. The speed of heat 
propagation diverges in the limit of vanishing r and be¬ 
comes infinite for the system described by Fourier’s law. 

It is a common understanding, that Eq. HD, being 
treated in the context of CIT, leads to violation of the 
second law of thermodynamics 0. It can be shown that 
the classical specific entropy oscillates with the course of 
time, if the system is described by Eq. Jl]). In particular, 
classical specific entropy can decrease with time, which 
contradicts to the requirement of the entropy produc¬ 
tion being always positive. To overcome this issue, there 
have been developed a theory known under the name ” ex¬ 
tended irreversible thermodynamics” (EIT) E3- The key 
difference of EIT with respect to CIT is the relation to the 
local equilibrium hypothesis. It is assumed in CIT that 
the thermodynamic relations formulated for equilibrium 
macroscopic systems are also valid in non-equilibrium, if 
applied locally in space and time. In other words, every 
small sub-volume of a system for a given moment of time 
is assumed to reach it’s equilibrium, such that there exist 
local equilibrium. Mathematically this is formulated such 
that the specific entropy is a function of the local values 
of the state variables (e.g. the local densities) only. In 
contrast, EIT assumes that in non-equilibrium the spe¬ 
cific entropy may in addition depend on the fluxes as 
independent variables. Therefore, the local thermody¬ 
namic relations in non-equilibrium are not the same as 
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in equilibrium, and the system is said to be not in local 
equilibrium. The important consequence of considering 
the specific entropy to be dependent on both, the state 
variables and the fluxes, is that the entropy production 
is always positive. In addition, dependence of the fluxes 
leads to the equation of the MCV type 0- Thus, consis¬ 
tency with the second law of thermodynamics is restored. 

In the earlier work 0 we have shown that the second 
law of thermodynamics can be viewed not as an indepen¬ 
dent law, but as a consequence of a variational procedure, 
which minimizes a certain thermodynamic ” action”. The 
important part of that procedure was the assumption of 
local equilibrium. It was possible to derive the classical 
force-flux relations, one of which had the form of Fourier’s 
law of heat conduction. In view of that, it is interesting 
to realize, whether a similar procedure can result in a 
force-flux relation of the MCV type 0 . Furthermore, it 
is interesting to understand whether local equilibrium is 
sufficient to describe the systems, which are governed by 
Eq. 0 . In this paper we will address these questions. We 
will show that the property of local equilibrium is, in fact, 
sufficient for the system to exhibit a MCV-like behavior. 
We will also address the issue of contradiction between 
the local equilibrium hypothesis and the second law of 
thermodynamics. In particular, we will show, that the 
formulation of the second law of thermodynamics, being 
considered as a consequence of the variational procedure, 
should be modified. 

Variational methods in irreversible thermodynamics 
have mostly been considered for restricted sets of con¬ 
ditions 00- In particular, they have mostly been 
formulated either for the constant values of the trans¬ 
port coefficients or for stationary conditions. The cel- 
ebratedprinciple of minimum total entropy production 
[p. 45] [ll, formulated by Prigogine, states that in sta¬ 
tionary states the evolution of the system with constant 
transport coefficients is such that the total entropy pro¬ 
duction is minimum. There have been other variational 
formulations of irreversible thermodynamics with vari¬ 
ous restrictions 0 , 0 |. Not much work has been done 
for the systems, which possess thermodynamic inertia. 
The systems of heat wave propagation, which are p urely 
reversible, allow a standard variational formulation [lj] . 
The MCV-like systems can be viewed as waves with dissi¬ 
pation, and the dissipation can be included by introduc¬ 
ing an exponentially decaying factor in the Lagrangian 
|15| . Such variational description can not be considered 
general enough, since the Lagrangian explicitly contains 
a particular solution. In addition, it is not symmetric 
with respect to time reversal, and therefore is not com¬ 
patible with microscopic time reversibility. 

In our earlier work [llj| we formulated a variational pro¬ 
cedure, which allows one to derive the force-flux relations 
and the second law of thermodynamics for the systems 
without thermodynamic inertia. We suggested the La¬ 
grangian which is symmetric with respect to time rever¬ 
sal and does not contain explicitly nether the solutions, 
nor the evolution equations. The variational principle is 


formulated for the so-called extended system, which con¬ 
sists of the normal system, the so-called mirror-image 
system, and the equilibrium system. Evolution of the 
normal system is characterized by the normal second law 
of thermodynamics, which states that during the evolu¬ 
tion the entropy production is always positive, cr > 0. In 
contrast, evolution of the mirror-image system is char¬ 
acterized by the negative entropy production, cr* < 0, 
which can be viewed as the mirror-image second law of 
thermodynamics. Finally, the equilibrium system is re¬ 
quired as a reference system for the other two. The de¬ 
scription of the irreversible evolution of such system is 
completely equivalent to the one, given by the standard 
CIT [lj. In this paper we extend this analysis to the 
systems with non-zero thermodynamic inertia. 

The paper is organized as follows. In Sec. [II] we briefly 
mention the important steps of the variational procedure 
for the system governed by CIT. In this way we also in¬ 
troduce the notation and the important principles. In 
Sec. HU we extend the variational procedure to the sys¬ 
tems with thermodynamic inertia. We present the La¬ 
grangian, which leads to MCV-like equations of evolu¬ 
tion. We also present a conserved quantity, which is 
an analogue of the Hamiltonian in classical mechanics. 
Conservation of this Hamiltonian during evolution is im¬ 
portant for understanding the entropy transformations in 
the system. In Sec. |TV]we consider the specific case of the 
MCV-like system, the thermal wave propagation. Such 
system has no dissipation but still has thermodynamic 
inertia. It is therefore convenient to view such system as 
a reference system for our analysis. In Sec. [V] we intro¬ 
duce the coordinates, in which evolution of the system 
has simplest form and is easy to analyze. We discuss the 
entropy transformation in the system with both non-zero 
thermodynamic inertia and dissipation in Sec. lVIl In par¬ 
ticular, we give a new formulation of the second law of 
thermodynamics, which does not contradict to local equi¬ 
librium. Finally, in Sec. IVIII we discuss the similarities 
between the presented theory and the known and devel¬ 
oped fields of physics, in particular, classical mechanics 
and electrical engineering. 

II. CLASSICAL IRREVERSIBLE 
THERMODYNAMICS 

Irreversible evolution of the system, which is described 
by the Fourier-like law between the thermodynamic force 
and the thermodynamic flux is derived from the principle 
of stationary action, with the Lagrangian given by 

L C1T = -~t{v eq ) V<p • Vp* - i (<pp* - <p*p) (2) 

Here and further the quantities with the tilde ~ will 
denote the properties related to the extended system, 
i.e. to the all three systems together, the normal one, 
the mirror-image one and the equilibrium one. The La¬ 
grangian of the extended system depends on the poten¬ 
tial ip and the material density p, which are related by 
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an equation of state p = p EOS (p)- In the case of heat 
conduction, the potential is 1/T, where T is the temper¬ 
ature, while the density is the internal energy density. 
In the case of diffusion the potential is —p/T, where p 
is the chemical potential, while the density is the den¬ 
sity of the component. In the case of electric conduction 
the potential is —< f>/T , where </> is the electric potential, 
while the density is the charge density. The important 
property of the Lagrangian © is that it is symmetric in 
time. The actual evolution of the system, which makes 
the irreversible action f L dr dt to reach its extremum, 
is described by the solution of the Euler-Lagrange equa¬ 
tions. 

The Euler-Lagrange equations represent the pair of 
the force-flux relations, in the normal system and in the 
mirror-image system: 


Lagrangian is performed. Because of this (. and x may 
have any dependence on the equilibrium state functions 
(e.g. the temperature), and do not have to be constants. 

We have also shown that the above procedure allows 
one to derive the second law of thermodynamics. The 
variational principle implies that in the course of evolu¬ 
tion there exists a quantity W, which is the analogue of 
the Hamiltonian density in classical mechanics: its inte¬ 
gral over the volume of the system is conserved with the 
time 

j /» 

— / nw = o (6) 

dt Jy 

The Hamiltonian density for the Lagrangian © has the 
form 

Wcit = £{p eq ) Vp • Vtp* (7) 


J= £Vp 
J* = -£Vp* 


( 3 ) 


The coefficient £ is the ordinary phenomenological coef¬ 
ficient jlj. It is positive. It is an even function with 
respect to time reversal. And it is independent of ei¬ 
ther p or p*. According to the Green-Kubo relations, it 
is equal to the equilibrium time-correlation function of 
corresponding microscopic fluxes. For not far from equi¬ 
librium perturbations we also have 

e(v )« ^ eq )« %*) ( 4 ) 

The actual evolution of the system obeys the prop¬ 
erty of local equilibrium. This, in particular, means 
that non-equilibrium equation of state has the same form 
as in equilibrium, p(r,t) = p EOS (p{r,t)) and p*(r,t) = 
/) EOS (p*(r,f)), where the function p EOS is defined by the 
equilibrium relation p eq = p EOS {p eq )- Because of this, 
we can write that 

dp*{r,t) = dp eq = dp{r,t) _ 

dp*(r,t) dp eq dp(r ,t) ^ 

where for further convenience we introduced a positive 
quantity x-, which we will call a capacity. Indeed, in the 
case of heat conduction the potential ip = 1 /T and the 
material density p = u is the internal energy density. 
The derivative du/d(l/T) = — T 2 C , where C is the heat 
capacity. The heat capacity is positive, which makes x 
to be positive and proportional to the ordinary heat ca¬ 
pacity. A similar identifications can be made in the case 
of diffusion and electric conduction. 

We should emphasize, that we require the property 
of local equilibrium to hold only for the actual evolu¬ 
tion trajectory, i.e. the one, which satisfies the Euler- 
Lagrange equations. This means that Eq. © and Eq. © 
are correct only for the solution of the Euler-Lagrange 
equations, but not for any other evolution trajectory. In 
other words, the assignments © and a should not be 
made directly in the Lagrangian, but only in the resulting 
Euler-Lagrange equations, i.e. after the variation of the 


which implies that the entropy productions in the nor¬ 
mal and the mirror-image systems have opposite signs. 
Indeed, their product is negative, era* = -tV E i T , and 
the entropy production of the normal system is always 
positive, a > 0. 


III. MAXWELL-CATTANEO-VERNOTTE 
EQUATION 


Let us consider the following Lagrangian 

Lmcv = m{p eq ) pp* - £(p eq ) Vp ■ \7p* - i (pp* - p*p) 

( 8 ) 

It differs from L CIT by an additional term, which is pro¬ 
portional to the product of the time rates of change of the 
material densities of in the normal and the mirror-image 
systems. As we will see further, the proportionality co¬ 
efficient m is a mere of the thermodynamic inertia of the 
system. Just like the coefficient £, m is positive and is 
an even function with respect to time reversal, and also 
is independent of either p or ip*. The Lagrangian © is 
symmetric in time, just like the Lagrangian (©. 

Just like for the classical irreversible thermodynam¬ 
ics, we assume that the actual evolution of the system 
obeys the property of local equilibrium. This means that 
Eq. © and Eq. © are still true for the trajectory, which 
makes the irreversible action J L MCV dr dt to be extremal. 
In addition, we require for the extremal trajectory not far 
from equilibrium that the inertial function m satisfies 


m(p) ~ m{<p eq ) ~ m(p*) (9) 

The Lagrangian © depends on the first-order spatial 
and temporal derivatives of the state functions. The 
Euler-Lagrange equations have therefore the following 
form 


dLyiCV 

dp 


9L MC y 

dp* 


d dL MCV 
dt dp 

d dL MCV 
dt dp* 


dL M cv 

dS7p 

9LmCV 

dWp* 


= 0 

= 0 


( 10 ) 
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Evaluating the partial derivatives and taking into ac¬ 
count Eq. ©, Eq. © and Eq. I© we obtain the evolution 
equations in the normal and the mirror-image systems 


rp + p + V • (£ V<^) =0 

r/S*-/>*+V-(£V<p*) =0 


( 11 ) 


course of evolution there exists a Hamiltonian density 
W CIT , integral of which over the volume of the system 
is conserved with time, so that Eq. © holds for W C1T . 
The Hamiltonian density is W = — L + ip(dL/d<p) + 
(fi*(dL/dip*). In the case of Lagrangian ®, this results 
in the following expression 


where we introduced 


T{(p eq ) = xm ( 12 ) 


The parameter r is positive and has dimensionality of 
time. The terms rp and rp* represent delayed responses 
of the material densities to the perturbations Vtp and 
\7ip*. Thus, the parameter r may be identified with the 
delay time of the system’s response, which is a conse¬ 
quence of non-zero thermodynamic inertia. The delay 
time is proportional to the inertial parameter m and char¬ 
acterizes thermodynamic inertia of the system. In partic¬ 
ular, when the inertial parameter m = 0, the system has 
no inertia, and the evolution equations become equiva¬ 
lent to the force-flux relations of the classical irreversible 
thermodynamics. 

We next combine Eq. m with the balance equations 
for the normal and the mirror image convection-free sys¬ 
tems. There the sum of the divergence of the material 
flux and the time rate of change of the material density 
is equal to zero for both, the normal system, V J + p = 0, 
and the mirror-image system, VJ* + p* = 0. It follows 
therefore that 


r j + J — l\7ip = 0 
Tj*-J*-£V<p* =0 


(13) 


Eq. (fT3l) are the force-flux relations for an irreversible 
system with non-zero thermodynamic inertia. For heat 
conduction the potential ip = 1/T and it is easy to see 
that the first of Eq. m is the MCV equation m with 
the thermal conductivity related to the phenomenological 
coefficient in a standard way, A = £/T 2 . 

It follows therefore that evolution of the inertial sys¬ 
tem, which is described by the MCV-like equation, can 
be derived from the variational principle. The important 
step in the derivation is the property of local equilibrium, 
which allows one to apply Eq. © , Eq. © and Eq. © for 
the extremal evolution trajectory. The variational pro¬ 
cedure for the inertial system is the same as the one for 
the inertia-less system, introduced in El and outlined in 
the previous section. The difference between the evolu¬ 
tion of the inertial and inertial-less system is due to the 
additional term mpp* in the Lagrangian L Mcv compared 
to the Lagrangian L CIT . We shall call this term the ki¬ 
netic entropy production. When the inertial parameter 
m = 0, the Lagrangian L MCV becomes equivalent to the 
Lagrangian i CIT , and the evolution of the system follows 
the classical Fourier-like equation of heat conduction. 

As in the case of inertia-less system, the variational 
procedure for the inertial system implies that in the 


W MCV =%> e9 )Vv?-V<p* + m(<p eq )pp* (14) 

The Hamiltonian density for the inertial system differs 
from the Hamiltonian density for the inertia-less system 
by the kinetic entropy production. The conserved quan¬ 
tity in the inertial system, the volume integral over the 
Hamiltonian density, contains therefore a kinetic contri¬ 
bution. 


IV. RESISTLESS EVOLUTION 


Depending on the values of the parameters r and £, 
the system may exhibit specific behavior. In particular, 
when t is very large, the second term in Eq. m may 
be neglected, and we obtain two wave equations. Wave 
propagation is a conservative process, so the limit of large 
r corresponds to a system without dissipation, i.e. re¬ 
sistless system. However, the parameter r characterizes 
delayed response of the system to external perturbation, 
which seems to be unrelated to the dissipative properties 
of the system and it might be not convenient to use r as 
a mere of dissipation in the system. We formalize this 
observation, we shall view the resistless system from a 
different perspective. 

In particular, let us introduce the resistivity r(tp eq ) = 
l/£. Furthermore, let us introduce the heat wave factor 
■w(ip eq ) and can the heat wave speed c(ip eq ): 


w 


2 


t_ _ 1 i 
m x T 


c 


2 


2_t_ 

X 3 T 


(15) 


As wave propagation is a conservative process, its speed 
may not depend on the irreversible characteristics, in par¬ 
ticular the phenomenological transfer coefficient l. We 
may expect therefore m to be proportional to i. Fur¬ 
thermore, we may expect that it is w (or c) that reflects 
the ’’true” inertial characteristics of the system, while m 
is a combination of the irreversible characteristic £ and 
the inertial characteristic w (or c). 

With the new quantities the Lagrangian © can be 
written as 


Amcv = l{v eq ) 




~\ r {V eq ) (W* - P*P) 


( 16 ) 
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It leads to the same irreversible evolution equations m 
and results in the same MCV-type equation (fldl) . The 
wave equation is obtained when r = 0. Indeed, in this 
case the Lagrangian becomes 


£ w = £(f eq ) 


■ pp* -S7p-S7p* 


(17) 


If the coefficients are constant, the minimizing the La¬ 
grangian L w results in the following wave equations, 
which have the same form for both, the normal system 
and the mirror-image system: 


ip = c 2 \/ 2 p 

(p* = c 2 vv 


(18) 


In other words, the resistless evolution corresponds to the 
wave propagation. We can see, that considering a system 
with a large relaxation time r is equivalent to considering 
a system with a small resistivity r. In particular, when 
r = 0, t is infinitely large. 

The Hamiltonian density yV w for the resistless system 
has the same form as for an ordinary MCV-like system, 
Eq. (HD: 


W w = m{f eq )pp*+£(ip eq )Vf Vf* (19) 

This means, in particular, that if a resistless system has 
the same initial and boundary conditions as the corre¬ 
sponding MCV-like system, then the conserved Hamilto¬ 
nian has the same value in both systems. 

It is interesting to observe, that the same wave equa¬ 
tions can be obtained from the Lagrangians, which in¬ 
volve the state variables from only one system, either the 
normal one, or the mirror-image one. Indeed, consider 
the following single-system Lagraninans 

^l a -|v*l a ' 

( 20 ) 

-4lP*l 2 -|W*| 2 

w z 

The Lagrangian L w describes evolution of the normal 
system only. Similarly, the Lagrangian L ((, describes evo¬ 
lution of the mirror-image system only. Performing the 
variational procedure for each of this Lagrangians inde¬ 
pendently results in the wave equation. In case of the 
constant coefficients they have the form of Eq. (fTBl) . As 
expected, they have the same form in both, the normal 
system and the mirror-image system. One can conclude, 
that evolution of the extended resistless system is equiv¬ 
alent to the combined evolution of the normal resistless 
system and the mirror-image resistless system. In partic¬ 
ular, the Lagrangian L w + L^ describe the same resistless 
system as the Lagrangian L w . 

The single-system Hamiltonian densities for the nor¬ 
mal and the mirror image systems have the following 


£{<p eq ) 


£{p eq ) 


form 


yw 


m(p eq ) 

2 

m(p eq ) 


IpI 2 


\ p *\ 2 + 


£{p eq ) 

2 

£{p eq ) 


| V ^| 5 


|VvV 


( 21 ) 


As W w , they satisfy the conservation equation It 
follows, in particular, that W w + = W w . 

We see that evolution of the extended system, which 
is resistless, can be considered as a joint evolution of the 
normal and the mirror-image system, which evolve in¬ 
dependently. In other words, resistless system is decou¬ 
pled into the normal one and the mirror-image one. It 
follows therefore that the resistivity r is a mere of cou¬ 
pling between the normal and the mirror-image system. 
When r = 0 these two systems are completely decou¬ 
pled, and each of them represent a conservative system 
with wave propagation. Nonzero r introduces coupling 
between these system, which makes each of them irre¬ 
versible. In particular, the volume integrals of W w and 
WjJ) are no longer conserved. However, the extended sys¬ 
tem with nonzero resistivity still has a conserved quan¬ 
tity, the volume integral of W w . 


V. COMPLEX COORDINATES 


To understand the structure of, in particular, entropy 
transformation in the dissipative system with thermody¬ 
namic inertia, it is convenient to perform the analysis in 
complex space. We shall do that in this section. 


A. The generalized force 


A thermodynamic force, which drives the system away 
from equilibrium, is usually identified as A = S7p. In 
particular, for heat conduction it is V(l/T), for diffusion 
at constant temperature it is — V(///T), for electric con¬ 
duction at constant temperature it is —'V(<j>/T). Let us 
consider a complex force, which is a generalization of the 
standard thermodynamic force. In particular, we intro¬ 
duce 


3 = S7p + ie v — p (22) 

w 

where i = is the imaginary unit and is the unit 
vector in the direction of the gradient of ip. Together with 
5 we can consider the complex-conjugate force 5, the 
mirror-image conjugate force 5*, and their combination 
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S = V ip — ie v — p 
w 

E* =Vp*-ie v ,-p* (23) 

w 

2* = V<p*+ie lp * -p* 

w 

Using this notation, the Hamiltonian density for the 
MCV-like system can be written as 


with the coordinates |Vy>|, |Vy>*| and p, p*. Since the 
normal and the mirror-image systems are decoupled, evo¬ 
lution trajectory of the extended system is represented 
by two two-dimensional circles in the four-dimensional 
complex space, see Fig. 1(a) Each of this circle is 


a two-dimensional cross-section of the evolution trajec¬ 
tory of the extended system, which correspond to the 
pair of coordinates either (|Vv?|;p) or (|V(^*|; p*). The 
radii of these circles are constant and equal to |2| and 
|2*| respectively. Note, that for the resistless evolution 
Ww = W w + W^- This means that for the resistless 
evolution 2* = 2 and 5 = 2*. 



Furthermore for the resistless system, 


and 



while for the inertia-less system 


W c 


(24) 


(25) 


(26) 


(27) 


Let us mention first the inertia-less system. Evolution 
of such system was studied in [llj] and outlined in the 
second section of this paper. The generalized forces in 
this case is real, i.e. 2 = 2 and 2* =2*. The conserved 
quantity is 2 • 2*. 

Let us consider now the resistless system. As we saw 
in the previous section, evolution of the normal and the 
mirror-image system is decoupled. Individual evolution 
of the normal system can be represented by a trajectory 
in the two-dimensional complex plane with the coordi¬ 
nates |V</?| and p. Conservation of W w means that the 
absolute value |2| 2 of the generalized force is constant 
during evolution. Thus, the evolution trajectory of the 
normal system is a circle, and its evolution can be repre¬ 
sented as rotation of the vector 2 in the complex plane. 
The same arguments are valid for the mirror-image sys¬ 
tem as well. We have therefore 


P 

w| V tp\ 


= tan(cct) 


P* 

w|V<^*| 


= — tan(wi) 


(28) 


where w is the frequency and the sign of the right hand 
side is determined by the initial condition. On the other 
hand, evolution of the extended system, which contains 
the normal one and the mirror-image one, is represented 
by a trajectory in the four-dimensional complex space 



FIG. 1: Evolution trajectory in the complex space, (a) Re¬ 
sistless system, (b) Dissipative system 


Finally, let us consider the extended system with non¬ 
zero resistivity and non-zero inertia. The normal and 
the mirror-image systems are coupled, which means that 
both |2| and |2*| change. The two-dimensional cross- 
section of the evolution trajectory are now not circles, 
but spirals. The specific shape of the evolution trajectory 
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depends on the relationships between the parameters of 
the system and the boundary conditions. In particular, 
in the case of a relaxation process, Eq. m has a form of 
the equation for a damped oscillator with the damping 
ratio of the order ycr. When x cr ~ 1 th e normal sys¬ 
tem experiences decay with oscillations, while the mirror- 
image system expands with oscillations, see Fig. |l(b)| 
When x c r > 1 the corresponding relaxation and tight¬ 
ening processes happen without oscillations. In general, 
the vector of the generalized force changes such that the 
volume integral of the expression 5 ■ 3* + 3 • 5* remains 
constant. Also, 3* ^ 3 and S ^ 3*, as well as 3 ^ 3 and 
5* ^ 5*. This means that in general W M cv ^ W w + 
and, in particular, the value, the volume integral of which 
remains constant, is neither |S| 2 + |3*| 2 , nor 3 ■ 3*. 


The situations described above can also be illus¬ 
trated by the diagrams in Fig. [2] The resistless system 
(Fig. 2(c)) allows independent wave propagation in both 
the normal and the mirror-image system. The inertia-less 
system (Fig. |2(b)[ ) implies coupling between the normal 
system and the mirror-image one, as well as between their 
complex-conjugates. The full dissipative system with in¬ 
ertia (Fig. |2(a)j ) allows both coupling between the normal 
system and the mirror-image one, as well as wave prop¬ 
agation in either of them. The full system may not be 
decoupled in pairs like the resistless one or the inertia-less 
one. 





(b) 



FIG. 2: Interactions between the generalized forces in the (a) 
full system; (b) dissipative system; (c) resistless system. 


B. Impedance and force-flux relations 


Let us apply the Laplace transform to the MCV-like 
equation (flUl) . The Laplace transform of a function / is 
defined as /(k, oj) = exp {—wt) f dt. Here w = n + iou 

is a complex frequency, where k is the relaxation param¬ 
eter and u> is the wave frequency. The complex frequency 
characterizes the boundary (external) conditions, applied 
to the system (i.e. the temperature difference on the 
boundary, or the flux across the boundary). In particu¬ 
lar, if the system is subjected to a steady-state boundary 
conditions which are periodic in time with the period 
27 r/w, then vj = iu is purely imaginary. In contrast, if 
the system is subjected to a steady-state fixed boundary 
conditions, then vo = k is purely real. A transient or 
a relaxation process is characterized by a complex fre¬ 
quency, with neither k nor w are equal to zero. Laplace 
transform of Eq. (TTfll) has the following form: 

J(1 + r zu) + J(0) = 

(2D) 

J*(— 1 + ttu ) + J*(0) =£V(?* 


The complex impedance is defined as a ratio of the 
Laplace transform of the force and the Laplace transform 
of the flux, when the initial values of the flux is taken to 
be equal to zero. We have therefore 


Z(k,uj)= (1 + r k) r + iui t r 
Z*(k,lo)= (— 1 + Tn)r + iuiTr 


(30) 


The complex impedance for each of the system consists 
of the real ans imaginary parts. The real part determines 
the active resistance and characterizes dissipation in the 
system. The relaxation parameter for the normal system 
is equal to — r/r, which indicates exponential decay of 
a perturbation. In contrast, the relaxation parameter of 
the mirror-image system is equal to r/r, which indicates 
an exponential growth of a perturbation. The imaginary 
part of the complex impedance determines the reactive 
resistance. It does not lead to dissipation or accumu¬ 
lation, but characterizes conservative transformations of 
the entropy in a heat wave. The reactive resistance is the 
same for both, the normal and the mirror-image systems. 

Ignoring the initial condition, Eq. (l29l) can be written 
as 


V </? = Z(k , w ) J 
V £*= Z*(k,cu)J* 


(31) 


which is the force-flux relation for the irreversible system 
with thermodynamic inertia. 


VI. LOCAL EQUILIBRIUM AND THE 
ENTROPY PRODUCTION 

The derivations in this paper rely on the property of 
local equilibrium. The property of local equilibrium is 





























also crucial for classical irreversible thermodynamics. In 
particular, it is used when the Gibbs equation is formu¬ 
lated. This leads to the following expression for the local 
entropy production (which has the same form both for 
the normal and the mirror-image system): 


cr(r, t) = J ■ 
cr*(r, t) = J* ■ V<p* 


(32) 


If the system is inertia-less, the force-flux relations can 
be derived from the principle of stationary action [Tl] 
and are given by Eq. ©. Substituting these relations in 
Eq. (l32l) . leads to cr(r, t) > 0 and er*(r,f) > 0. The first 
of these relations is the statement of the second law of 
thermodynamics: the entropy production in the normal 
system is always positive, and reaches zero in equilibrium. 

If the system has non-zero thermodynamic inertia r, 
the force-flux relations can still be derived from the prin¬ 
ciple of stationary action, which is the subject of this 
paper. These relations are given by Eq. m instead 
of Eq. (|3j). Still, the derivation of Eq. (fl3l) implies the 
property of local equilibrium. Substituting Eq. m in 
Eq. (EEfl) results in the entropy production, which can¬ 
not, in general, be represented as either a positive or a 
negative quantity. This leads to apparent violation of 
the second law of thermodynamics for the system with 
thermodynamic inertia. 

There are two different conclusions which can be made 
on the basis of this observation: either the second law 
of thermodynamics should be reformulated, or the prop¬ 
erty of local equilibrium should be questioned. The latter 
conclusion is employed by extended irreversible thermo¬ 
dynamics [1]. In particular, the non-equilibrium Gibbs 
relation is modified such that the non-equilibrium local 
entropy depends not only on the state variables, but also 
on the irreversible fluxes. This results in the expression 
for the entropy production, which is always positive in 
non-equilibrium. 

The alternative conclusion, which requires revisiting 
the second law of thermodynamics, may sound unrea¬ 
sonable. Still, the analysis in this paper shows, that the 
MCV-like force-flux equations (1131) can be derived rig¬ 
orously within the assumption of local equilibrium. The 
property of local equilibrium is sufficient for deriving the 
MCV-like force-flux equations. In other words, there is 
no need to go beyond local equilibrium and assume addi¬ 
tional relationships. This suggests, that the second law of 
thermodynamics should be reformulated, such that evo¬ 
lution of an irreversible system with non-zero thermody¬ 
namic inertia would be in agreement with it. 

Eq. (l32l) describes the entropy production, which is lo¬ 
cal not only in space, but also in time, i.e. instantaneous. 
For the inertia-less system the instantaneous entropy pro¬ 
duction is a measure of instantaneous dissipation. How¬ 
ever, for the system with thermodynamic inertia, the sys¬ 
tem response to a perturbation is not instantaneous. In 
particular, the instantaneous flux depends not on the in¬ 
stantaneous force, but on the force at time r earlier. This 


leads to the observation, that the instantaneous entropy 
production is not very practical quantity to characterize 
dissipation in the system. 

To illustrate this, let us consider the resistless system, 
described in Sec. EY1 The evolution of the resistless sys¬ 
tem is a propagation of a heat wave. As every wave prop¬ 
agation, it is conservative. Still, this system obeys the 
property of local equilibrium, and the instantaneous en¬ 
tropy production for this system is given by Eq. (1321) . For 
the process of wave propagation, the instantaneous en¬ 
tropy production of the normal system oscillates around 
zero, taking both the positive and the negative values. 
Obviously, the instantaneous entropy production does 
not satisfy the second law of thermodynamics. Further¬ 
more, according to the standard understanding of the 
entropy production, the wave propagation should be a 
process with alternating energy dissipation and energy 
gain. However, it is neither of these two: wave propaga¬ 
tion is a conservative process. It follows therefore that 
instantaneous entropy production has no practical mean¬ 
ing for the resistless system. 

For the resistless system one can distinguish two kids 
of terms, which contribute to the Hamiltonian m- 
They have the dimensionality of the entropy produc¬ 
tion, so we will refer to them as such. Namely, one 
can speak of the kinetic entropy production m |p| 2 /2 and 
m 1/3*172 and the potential entropy production £ |V</3| 2 /2 
and £|V<p*| 2 /2. The kinetic entropy production exists 
due to the system’s inertia mneqO , while potential en¬ 
tropy production exists due to system’s inhomogeneity 
S/tpneqO. During the process of wave propagation the en¬ 
tropy is transferred from the kinetic term to the potential 
term and vise versa, keeping the Hamiltonians f W w dV 
and f dV constant. The instantaneous entropy pro¬ 
ductions (1321) also oscillate. This process is the same in 
both, the normal system and the mirror-image one. 

In contrast, the dissipative inertia-less system does not 
have the kinetic term. Furthermore, the potential en¬ 
tropy production has the form of £V(p ■ V(/?*, i.e. rep¬ 
resents coupling between inhomogeneities in the normal 
and the mirror-image systems. In other words, the po¬ 
tential entropy production is shared between the normal 
and the mirror-image systems. Still, the Hamiltonian of 
the extended system f W C i T dV is conserved. This makes 
the entropy to be transferred from the mirror-image sys¬ 
tem to the normal system. This process is irreversible, 
in contrast with the entropy transfer in the resistless sys¬ 
tem. 

The extended system with non-zero inertia and non¬ 
zero resistance exhibit both processes of the entropy pro¬ 
duction transfer: reversible transfer between the kinetic 
and the potential term, as well as irreversible transfer 
between the mirror-image and the normal system, - such 
that the Hamiltonian of the extended system J W MC v dV 
is conserved. It is clear, that the oscillatory part of the 
instantaneous entropy production is reversible and there¬ 
fore does not contribute to dissipative behavior. 

In order to separate wave propagation from dissipa- 
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tion, we should consider the entropy production, which 
is averaged over the wave period r w : 


1 r t+T " 

(er(r,f)) w =— / a(r,t')dt' (33) 

Tw J t 


and (cr*(r, t)) w is introduced in the same way. Note, that 
the wave period r w is different from the delay time r. 
It depends on the system’s size, as well as the speed 
of heat wave c. The statement of the second law of 
thermodynamics should then read, that the quantity 
which is positive in the course of evolution is the time- 
averaged entropy production of the normal system. Cor¬ 
respondingly, the time-averaged entropy production of 
the mirror-image system is negative: 


(cr(r,i)) w > 0 
(<r*(r,f)) w <0 


(34) 


Substituting Eq. (fT51) in Eq. (1551) we obtain the fol¬ 
lowing expressions for the time-averaged local entropy 
productions: 


(ct) w = r [A w (J 2 ) + (J 2 )] 
(ct*) w = r [A w (J* 2 ) - (J* 2 )] 


(35) 


where A w /(t) = f(t + r w ) — /(t). It is convenient to 
examine the validity of Eq. (1341) for partial solutions. 
For the resistless system r = 0 and both (cr) w = 0 
and (tx*) w = 0. For the steady heat wave caused 
by the oscillating external force or boundary conditions 
J oc exp(iwt), and therefore A W (J 2 ) = 0. This leads to 
(tr)w = r (J 2 ) > 0. Similarly, for the mirror-image sys¬ 
tem (<r*) w = —r (J* 2 ) < 0. For the stationary process J 
is independent of time, and therefore A W (J 2 ) = 0, which 
also leads to (cr) w > 0. Similarly, for the mirror-image 
system (tr*) w < 0. Finally, for the relaxation process, 
which is the result of overdamping in the oscillatory sys¬ 
tem, J oc exp(— nt) and (<r) w oc 1 — kt. In order for 
the time-averaged entropy production to be positive, we 
should have kt < 1. However, this is exactly the con¬ 
dition of overdamping, which is true for the considered 
relaxation process by definition. Similarly, J* oc exp (nt) 
and (cr)w oc —1 T kt < 0. 

We note that one should not confuse the instantaneous 
entropy production, which is used in this paper in the 
context of mesoscopic continuous description, with the 
instantaneous entropy production which is formulated 
with the help of microscopic distribution function. Cor¬ 
respondingly, one should not confuse time averaging over 
the period of the heat wave with time averaging over the 
time of molecular equilibration. The analysis of this pa¬ 
per is performed withing the assumption of local equilib¬ 
rium, which means that the system has already reached 
equilibrium locally, for a given position r in space and 
a given moment t in time. In particular, it is assumed 
that the time scale of molecular equilibration is much less 


than the period of the heat wave, so that molecular fluc¬ 
tuations have already been averaged at every moment t 
in time. 

Classical irreversible thermodynamics distinguishes 
between two types of entropy change [lj. The first one, 
d e S = 5Q/T , is the equilibrium entropy change due to 
reversible interaction with the environment. The second 
one, diS = / adV, is the production of entropy due 
to irreversible processes in the system. These transfor¬ 
mations can be illustrated by Fig. [51 where we can see 
the entropy transformations in the systems, consisting of 
two thermal baths and the environment. The baths are 
in equilibrium with environment, but not in equilibrium 
with each other. Such system is open, i.e. not conserva¬ 
tive, as the entropy is produced in the course of evolution. 



The above analysis shows, that the entropy transfor¬ 
mations in the irreversible system with non-zero ther¬ 
modynamic inertia are different from those mentioned 
above. In particular, we can distinguish between three 
types of entropy change in the normal system and three 
types of the entropy change in the mirror-image system: 


dS = d e S + diS + d m S 
dS* = d e S* + d>S* + d m S* 


(36) 


The terms d e S and diS in the entropy change of the nor¬ 
mal system are the standard terms of classical irreversible 
thermodynamics. The third term, d m S 1 is present due 
to non-zero thermodynamic inertia of the system. This 
term is reversible and accounts for the entropy transfer in 
the heat wave. Similarly, we can speak of the ’’classical” 
terms d e S* and diS* in the mirror-image system. d e S* 
is responsible for reversible exchange with environment, 
while diS* represent irreversible entropy destruction. In 
addition, there exists the third term d m S *, which has 
the similar nature as d m S. In particular, it is reversible 
and responsible for the entropy transfer in the heat wave 
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in the mirror-image system. These entropy transforma¬ 
tions can be illustrated in Fig. [4j The extended system 
is closed. Indeed, the entropy transfer caused by the 
thermodynamic inertia is reversible and confined within 
each of the two subsystems, either the normal one or the 
mirror-image one. The irreversible entropy transfer hap¬ 
pens from the mirror-image system to the normal system, 
i.e. is not produced or destroyed anywhere beyond the 
extended system itself. 



FIG. 4: Entropy transfers in the extended system. 


VII. DISCUSSION 


In this paper we considered evolution of the irreversible 
system with non-zero thermodynamic inertia. Behavior 
of such system has many analogies to the behavior of the 
systems, studied within well-established areas of physics, 
where inertia is important. In particular, these areas are 
classical mechanics and electrical engineering. 


A. Classical mechanics 


The analogy with classical mechanics can be evident 
from the terminology used in the paper. Thus, for the 
resistless system the inertial entropy production mp 2 /2 
is the direct analogy of the kinetic energy mi /2 of a body 
with mass m and position x , while the potential entropy 
production £|V<^| 2 /2 is the direct analogy of the poten¬ 
tial energy k( Ax) 2 /2 of the body at the position Ax away 
from equilibrium with the force constant k. The potential 
energy of a body in classical mechanics is a function of 
its position only, and is independent of the velocity. Sim¬ 
ilarly, the potential entropy production is a function of 
the thermodynamic potential ip only, and is independent 
of the its ’’speed” p or the flux J. Thus, the property of 
local equilibrium in irreversible thermodynamics is nat¬ 
urally satisfied in classical mechanics. The sum of the 
potential and the kinetic energy of an isolated mechani¬ 
cal system is conserved. Similarly, the sum of the poten¬ 
tial and kinetic entropy production of a resistless ther¬ 
modynamic system is conserved, which is illustrated, in 
particular, in Fig. 1(a) Furthermore, the analogy of the 


heat wave factor w in classical mechanics is the eigenfre- 
quency uj. The wave factor does not have dimensionality 
of the inverse second, because we consider here a continu¬ 
ous system, which implies wave propagation, rather than 
oscillation of a single body. The typical characteristics of 
such system is the speed of wave propagation, c, which 
in our case is the speed of heat wave. 

The thermodynamic system with non-zero resistivity is 
analogous to the mechanical system with non-zero fric¬ 
tion. Conventional Hamiltonian formalism in mechanics 
is not applicable for such systems. However, it is still 
possible to account for friction by introducing the mirror- 
image system for a mechanical system as well [l 4| . The 
mirror-image system gains the energy dissipated from the 
normal system. Similarly, in irreversible thermodynam¬ 
ics, the mirror-image system gains the entropy, which is 
produced in the normal system, Indeed, the variational 
approach, which is applied in this paper, implies exis¬ 
tence of a quantity, which has a minimum for the actual 
evolution, the action. For a thermodynamic system it has 
dimensionality of the entropy generation, and it is nat¬ 
ural to understand the action as a kind of entropy gen¬ 
eration of the extended system (which consists of both, 
the normal one and the mirror-image one). The varia¬ 
tional formalism implies the property, which can be for¬ 
mulated as following: the total entropy generation of the 
extended system during evolution reaches its minimum. 
For the resistless system this leads to the conservation of 
the normal entropy, while for the dissipative system this 
leads to the entropy transfer from the mirror-image sys¬ 
tem to the normal system. The important consequence 
of this is that the normal and the mirror-image system 
should be considered together and some of their prop¬ 
erties may be shared within the extended system. The 
obvious example of the shared property is the Hamilto¬ 
nian. Furthermore, we can speak of the potential entropy 
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generation and the kinetic entropy generation 

mpp*, sum of which is equal to the Hamiltonian. Thus, 
these terms may also be viewed as analogies to the po¬ 
tential and kinetic energy in classical mechanics. 

The analogy between the mechanical and the thermo¬ 
dynamic quantities leads to analogy between the behav¬ 
ior of the mechanical and thermodynamic systems. A 
body with zero mass, which moves in media with friction, 
will eventually stop. If the body is attached to a spring 
and perturbed from its equilibrium position, it will stop 
in the spring equilibrium position, which corresponds to 
the minimum potential energy of the spring. This is a 
consequence of the so-called principle of minimum po¬ 
tential energy. Similarly, an irreversible thermodynamic 
system perturbed from equilibrium will experience re¬ 
laxation. Eventually it will reach the equilibrium state, 
which corresponds to the state with maximum entropy. 
This process is formulated as the second law of thermo¬ 
dynamics. If follows therefore, that the second law of 
thermodynamics is the thermodynamic analogy of the 
mechanical principle of minimum potential energy. 

The principle of minimum potential energy is violated 
for systems with non-zero mass. Indeed, if the body 
has non-zero kinetic energy, its potential energy may in¬ 
crease. For the system without friction this leads to the 
decrease of the kinetic energy, such that the total en¬ 
ergy remains the same. For the system with friction, the 
kinetic energy still can be transformed to the potential 
energy. The total energy, however, is not conserved, and 
the total energy averaged over the period of oscillations, 
decreases with time. Both systems, with and without 
friction, can be described by the variational formalism 
such that the energy is transferred between the kinetic 
and the potential form as well as between the normal 
and the mirror-image systems. Clearly, the principle of 
minimum potential energy is not directly applicable for 
the mechanical system with non-zero mass. Instead, one 
should use the principle of stationary action, either for a 
frictionless system or for a system with friction. 

Similarly, the second law of thermodynamics is not di¬ 
rectly applicable for a thermodynamic system with non¬ 
zero inertia. Instead, one should use the principle of sta¬ 
tionary action, which has been formulated in this paper. 
This principle leads to conservation of the total entropy 
generation in the extended system. 


B. Electrical engineering 


The analogy with the electrical engineering is evident 
from the fact, that MVC equation has the form of the 
telegrapher’s equation for an electrical transmission line. 
In particular, the pair of equation for the voltage V and 
the current / on the line, consisting the resistance R , the 


inductance L and the capacitance C is 
W = - Li-RI 

V/= -CV 


(37) 


The first of these equations is the analogy of MCV equa¬ 
tion m with 1/R being the thermal conductivity A and 
L/R being the inertial parameter r, while the second 
of these equations is the analogy of the energy balance 
equation, with C being the heat capacity. Furthermore, 
the analogy of the current / is the flux J, while the anal¬ 
ogy of the voltage V is the temperature T. The current 
propagation in such systems has been extensively stud¬ 
ied elsewhere. In particular, the inertial properties of the 
circuit are associated with a non-zero inductance L. Be¬ 
cause of it the sinusoidal voltage and sinusoidal current 
have different phase: the AC voltage leads the AC cur¬ 
rent. If there is no external electromotive force, the cur¬ 
rent decays with time. The active resistance is the source 
of energy dissipation, while the presence of inductor and 
capacitor introduce the reactive resistance. It does not 
lead to dissipation of energy, but allow the elements to 
gain and release the energy in a reversible way. 

The analysis above is directly applicable to such sys¬ 
tem. Since the temperature of the electrical circuit is 
considered to be constant and equal to the temperature 
of the environment, energy dissipation is proportional to 
the entropy dissipation. All the above relations and con¬ 
clusions remain valid if we simply replace the word ’’en¬ 
tropy” with ’’energy”. In particular, we should be able 
to distinguish between the kinetic and potential energy 
of the circuit. Furthermore, the impedance introduced 
in Sec. El is the electrical impedance, and the complex 
force-flux relations are the relations between the complex 
current and the complex voltage. 

The power of the electric circuit is equal to the prod¬ 
uct of the current and the voltage. For the inertia-less 
circuit all power is positive and active, i.e. is released on 
the resistor. For the circuit with non-zero inertia the in¬ 
stantaneous power does not have much physical meaning. 
Because the current and the voltage are not in phase, the 
instantaneous power does not have definitive sign and os¬ 
cillates. One can distinguish between the active power, 
which is released on the resistor, and the reactive power, 
which characterizes the oscillations of the electric cur¬ 
rent. To discard the power oscillations during one pe¬ 
riod of AC, it is useful to consider the power, which is 
averaged over the period of AC. The averaged power in¬ 
dicates the losses on the resistor only, and therefore is 
always positive. 

Comparing the behavior of the electric circuit and the 
inertial irreversible system, we can observe that the elec¬ 
tric power is the analogy of the entropy production. For 
DC current or for an inertia-less circuit the instantaneous 
electric power is always positive, which is the statement 
of the second law of thermodynamics. In contrast, for 
AC current in the circuit with inertia it is only the time- 
averaged power which is positive, while the instantaneous 
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power does not have definitive sign. The same applies to 
the entropy production in the electric circuit and in the 
irreversible system with inertia in general. 

C. Concluding remarks 

Let us summarize the findings presented in this paper. 
We have shown that irreversible evolution of a thermody¬ 
namic system can be derived from a variational principle 
of stationary action. In order to formulate this principle 
we need to consider the extended system, which consists 
of the normal system and the mirror-image system. The 
action, which is being minimized has the dimensionality 
of entropy and can be viewed as the entropy generation 
of the extended system. Thus, the principle of station¬ 
ary action states that evolution of the extended system 
minimizes the generated entropy of the extended system. 

The Lagrangian introduced in this paper describes gen¬ 
eral evolution of an irreversible system. It can also de¬ 
scribe the specific cases: resistless evolution of the sys¬ 
tem with non-zero inertia, which leads to wave propaga¬ 
tion; as well as dissipative evolution of inertia-less sys¬ 


tem, which leads to relaxation. The evolution equations 
or the force-flux relations are derived from the variational 
procedure. They represent the Euler-Lagrange equations 
for the corresponding Lagrangian. The important part 
of the analysis is the property of local equilibrium, which 
states that the non-equilibrium relation between the ma¬ 
terial density and the potential has the same form as in 
equilibrium. It appears, that the property of local equi¬ 
librium is sufficient for derivation of the experimentally 
verified evolution equations. Thus, no additional rela¬ 
tions which characterize non-equilibrium evolution are 
needed. 

Earlier it has been shown that dissipative evolution of 
inertia-less system results in the second law of thermody¬ 
namics. It states that instantaneous entropy production 
of the non-equilibrium normal system is always positive. 
Here we have shown that this formulation is no longer 
valid for the dissipative system with non-zero inertia. 
However, one can modify the statement of the second 
law of thermodynamics. In particular, it should read: 
for the irreversible evolution of the system the averaged 
over the period of the heat wave entropy production is 
always positive. 


[1] S. R. de Groot and P. Mazur. Non-Equilibrium thermo¬ 
dynamics. Dover, New York, 1984. 

[2] Peter Van and Tamas Fiilop. Universality in heat con¬ 
duction theory: weakly nonlocal thermodynamics. Ann. 
Phys. (Berlin), 524:470478, 2012. 

[3] Kumar K. Tamma and Xiangmin Zhou. Macroscale and 
microscale thermal transport and thermo-mechanical in¬ 
teractions: some noteworthy perspectives. Journal of 
Thermal Stresses, 21:405-449, 1998. 

[4] D. Y. Tzou. Macro-to-Microscale Heat Transfer. The 
Lagging Behaviour. Taylor & Francis, New York, 1997. 

[5] D. D. Joseph and Luigi Preziosi. Heat waves. Rev. of 
Mod. Phys., 61:41-73, 1989. 

[6] J. C. Maxwell. Philos. Trans. R. Soc. Lond., 157:49, 1867. 

[7] C. Cattaneo. Atti Sem. Mat. Fis. Univ.Modena, 3:83, 
1948. 

[8] M. P. Vernotte. C.R. Hebd. Seances Acad. Sci., 246:3154, 
1958. 

[9] David Jou, Jose Casas-Vazquez, and Georgy Lebon. Ex¬ 
tended Irreversible Thermodynamics. Springer, New York 
Dordrecht Heidelberg London, 4rd english edition, 2010. 


[10] G. Lebon and D. Jou. Early history of extended ir¬ 
reversible thermodynamics (1953-1983): An exploration 
beyond local equilibrium and classical transport theory. 
Eur. Phys. J. H, 40:205240, 2015. 

[11] K. S. Glavatskiy. Lagrangian formulation of irreversible 
thermodynamics and the second law of thermodynamics. 
J. Chem. Phys., 142:204106, 2015. 

[12] J. Schechter. The Variational Method in Engineering. 
McGraw-Hill, New York, 1967. 

[13] Tamas Matolsci, P. Van, and J. Verhas. Fundamental 
problems of variational principles: objectivity, symme¬ 
tries and construction, volume Chapter 3 of Variational 
and Extremum Principles in Macroscopic Systems. Else¬ 
vier, Amsterdam, 2005. 

[14] Philip M. Morse and Herman Feshbach. Methods of The¬ 
oretical Physics. McGrow-Hill, New York, 1953. 

[15] S. Sieniutycz. Variational approach to extended irre¬ 
versible thermodynamics of heat and mass transfer. J. 
Non.-Equilib. Thermodyn., 9:61-70, 1984. 



